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Abstract 

The von Neumann evolution equation for density matrix and the Moyal equation for the Wigner 
function are mapped onto evolution equation for optical tomogram of quantum state. The connec- 
tion with known evolution equation for symplectic tomogram of the quantum state is clarified. The 
stationary states corresponding to quantum energy levels are associated with the probability rep- 
resentation of the von Neumann and Moyal equations written for the optical tomograms. Classical 
Liouville equation for optical tomogram is obtained. Example of parametric oscillator is considered 



J> in detail. 

m 
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^ Introduction 



In [1], [2] new formulation of quantum mechanics called probability representation of quantum me- 
chanics (see recent reviews [3], [1]) was suggested. In the probability representation the quantum states 
are described by fair probability distributions called quantum tomograms or tomographic probability 
distributions. The tomograms contain complete information of the quantum state and they are con- 
nected with density operators of the states by means of invertable maps. There are several kinds of 
the tomograms related to the density operators by means of the different invertable maps. The idea of 
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introducing the probability representation of quantum mechanics was induced by existence of connection 
of the Wigner function W{q, p) of a system with continuous degrees of freedom hke position with its 
Radom transform [5j 6) called optical tomogram and found in [6], [?]. The optical tomogram and 

its connection with the Wigner function was applied in experiments on homodyne detection of photon 
quantum states [8] (see review article [S]), where the measured optical tomogram was considered as a 
technical tool to measure (reconstruct) the Wigner function of the photon state. 

In the approach based on reconstructing the Wigner function this quasidistribution function was 
considered as object identified with the quantum state and the optical tomogram was interpreted as 
intermediate step (not identified with the state) to obtain the Wigner function. In [10] the probability 
distribution called later [IT] symplectic tomogram was introduced. The symplectic tomogram M{X, fi, v) 
which is nonnegative probability distribution of homodine quadrature (of position) X depending on extra 
real parameters ^ and v was shown to be related with optical tomogram and this relation provided the 
possibility to reconstruct the Wigner function using the symplectic tomogram by means of Fourier 
transform of the tomogram. 

Since symplectic and optical tomograms are connected by invertable maps with the Wigner function 
and, consequently, with the state density operator it was suggested to identified the quantum states 
with the tomographic probability distributions [I], [2], [I2], [IS], [H]. It was shown that also in classical 
statistical mechanics the states can be described by the tomograms [12] , [13] connected with probability 
distributions f{q,p) on the system phase-space by the Radon transform. Thus, it was understood 
that the concept of state in both classical and quantum mechanics is identified with the tomographic 
probability distribution. The basic equation of quantum mechanics, like quantum evolution equation for 
density matrix [T3] and Moyal equation ^B] for the Wigner function [T7| were obtained in the probability 
representation for the symplectic tomograms [1], [2], [12], [13], [3]. Though the symplectic M(X, yU, z/) 
and optical tomogram w{X, 9) are connected by the relation w{X, 9) = M{X, cos 9, sin 9); the evolution 
equation known for symplectic tomogram was not transformed in explicit form into the evolution equation 
for the optical tomogram. Also the quantum equations for stationary states providing energy levels and 
written in the probability representation of quantum mechanics for symplectic tomograms (see, e.g. [3]) 
were not transformed till now into equation for the optical tomograms. On the other hand in experiments 
on homodyne detecting the photon states namely the optical tomogram is measured. In view of this 
one needs explicit form of the basic quantum equations written for the optical tomograms. Though, in 
principle, the equation written for symplectic tomogram implicitly contain information on the equation 
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written for the optical tomogram there is a technical difficulty to obtain the equation since the optical 
tomogram depends on two variables and the symplectic tomogram depends on three variables. One 
needs to take into account this difference to do the calculations. Thus the aim of our work is derive 
the evolution and energy level equations for the optical tomograms of quantum states for systems with 
several degrees of freedom and find the connection of these equations with equations for the symplectic 
tomograms. 

The paper is organized as follows. In Sec. 2 we review optical tomography of quantum states. In 
Sec. 3 and Sec. 4 the connection of the Moyal equation with the quantum evolution for optical tomogram 
and symplectic tomogram is studied and energy level equation for optical tomogram is presented. In 
Sec. 5 the classical Liouville equation for optical tomogram is derived. The tomograms of time dependent 
photon-added coherent states of parametric oscillator are found in Sec. 6. The conclusion and prospects 
are presented in Sec. 7. 



2 Symplectic and optical tomograms 

We review in this section the constructions of the optical and symplectic tomograms for system with 
one degree of freedom (or one- mode electromagnetic field). Given the state density operator p of the 
system. The symplectic tomogram M(X, /i, z/) is defined as 

M(X,/i, v) = (5(Xi -pq- up)) = — [ dk Tr{p e^'^^^-^"^"'^'^}. (1) 

2vr J 

Here q and p are position and momentum operators (photon quadratures), respectively, 1 is identity 
operator, and we assume that the taking trace and integration ower wave number k are permutative 
operations for the states under consideration. The optical tomogram w{X,6) is defined by analogous 
relation 

w{X,9) = {5{Xi~qcos9-psme)) = — [ dk Tr{p ei*^(^i-9=°««-P«i^e)}. (2) 

2tt J 

The phase < < 27r is called in homodyne detecting the photon states as local oscillator phase 
and this parameter can be controlled and varied in experiments on measuring the quantum states [H],[9]. 
Since the Dirac delta-function in Eq.([T]) is homogeneous function, i.e. S{Xy) = \X\~^6{y) the symplectic 
tomogram is the homogeneous function 

M(AX, A/i, Xu) = |A|-^M(X, p, v). (3) 
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From the property of delta function J S{y — a)dy = 1 it follows normalization condition for both sym- 
plectic tomogram and optical tomogram 



/ 



M(X,/i,z/) dX = 1 



(4) 



and 



jw{X,9)dX = l. (5) 
The homogenety condition ^ provides the relation of the optical and symplectic tomograms 

M{X, fx, v) = (/x^ + z/2)-i/2 w (x{fx^ + ^/^)"^/^ tan"^ 



V 

/i 

The Wigner function of state with density operator p is given by relation (we take h=X) 

u\ I u 



W{q,v) = J Trjp 
The Wigner function is normalized by the relation 



du. 



W{q,p) 



dq dp 

271 



(6) 



(7) 



(8) 



The symplectic tomogram is expressed in terms of the Wigner function in terms of its Radon transform 
depending on two parameters n and z/ 

dq dp 



M(X,/i,z/) 



W{q,p)5{X — /ig — up)- 



2% 



(9) 



The inverse transform is given by Fourier integral 

W{q,p) = — [ M(X,/i,z/)e*(^-^^-'^PMX d/x du. (10) 
2-71 J 

Analogously the optical tomogram is expressed in terms of Radon transform depending on the angle 6 

(11) 

Using (g one can to give the inverse Radon transform inserting the function M(X, /i, u) given by this 



wiX,9) = A /" H/(g,p)e^*^(^-9™«^-P«i"«)dg dp dA;. 
47r^ J 



relation into relation (|10|). Then one has 
W{q,p) 



1 



27r 



In this formula — oo < X, fi, u < oo. One can transform the integral using polar coordinates 



11 = r cos 6, v = r sin 6. 



(13) 



The symplectic tomogram is even function. It means that the optical tomogram has the property 



w{-X,e + 7i) =w{X,9). (14) 



Using this property the integral (12) can be presented in the form 

vr +00 +00 



W{q,p) = ^j de j y w;(X,^)|r/| e^''(^-''™^^-P^'^^Mr/ dX. (15) 







The optical and symplectic tomograms are probability distributions of random quadrature X. Thus one 
has two characteristic functions for these distributions 

XM{z,i^,iy) = J M{X,ij,u) e'^'^dX (16) 

and 

XUV,0) = J wiX,9)e^^''dX. (17) 
One has by using Fourier transforms the expressions 

M{X, /i, zy) = i- y xm{z, /i, ly) e-'^'^dz (18) 

and 

wiX, d) = ^J X^iv, 0) e-'^^'dv. (19) 
The characteristic functions are related 

Xw iv, 0) = XMiv, COS 9, sin 9) . (20) 

The highest momenta of the homodine quadrature {X"'){9) are determined by the characteristic function 
Xw{v,0), i.e. 

00 .„ 

xUv,0) = J2'^V''{X"m. (21) 

n=0 

3 Evolution equation for tomograms 

In this section we derive the evolution and energy level equations for optical tomograms. The von 
Neumann evolution equation for density operator p(t) reads 

dm 



dt 



+ i[H,p{t)] = ^- (22) 



Here the Hamiltonian H = + U{q) where U{q) is potential energy and the mass of the particle is 
taken as m = 1. In coordinate representation the von Neumann equation takes the form 



.dp{x, x', t) 



1 / 



dt 2 \dx^ (9x'2 

Using the relation of the density matrix and the Wigner function 



p(x, x', t) + (f/(x) - U{x')) p{x, x', t). 



p{x,x',t) = ^ I W 



X + x' 



,p,t] e 



^ip(x-x') 



dp 



one has the correspondence rules 

dp 



dW dp fl d \ 
-P^(^+25^j^; -P^[^--2d-v^'^- 



(23) 



(24) 



(25) 



which can be used to get the von Neumann equation (|23|) in the Moyal form 

i d 



dW{q,p,t) ^ dW{q,p,t) ^ 1 
dt dq i 



U\q 



2 dp 



c.c. 



W{q,p,t) = Q. 



(26) 



Using the relation (11) one can find correspondence rules for operators acting on the Wigner function 



and the optical tomogram. In fact 



d 



cos9-^w{X,e) = ^ j ikcosO e-''"^'°'''W{q,p,t) e^'^^^-P^'^'^^dk dq dp. 



(27) 



In view of the equality 



jfccos^ g-ifc<?cose 



d_ 
dq 



ikq cos 9 



and evaluating by parts the integral (27) over variables for functions W{q,p) — )■ for g — > ±00 we get 

dW{q,p) 



oos»A..(.Y.O) = ^ 



dq 



-5{X — q cos 6* — p sin 9)dq dp. 



It means that 



Q^W{q,p) M cos^— ti;(X, 



Analogously we have the correspondence rule 

d 



d 



^^Wiq,p)^sme—wiX,e). 
Applying the operator [d/dX)^^ which is defind by acting on plane wave as 



(28) 



(29) 



(30) 



_d_ 
dX 



-1 



ikX 



ikX 



ik 



(31) 



to (11) together with differentiation over the angle variable 9 and multiplication by sin 9 we get the 



equality 



89 [dJ^ 
We used identities 

and 



-1 



sin ^ j w{X,9)+X cos 9 w{X, ^) = ^ j Q W{q, p) e*fc(^-</cose-psine)^^ ^32) 



(q sin 9 — p cos 9) sin 9 + X cos ^ = g + (X — q cos ^ — p sin 9) cos ^ 



(5(X — q cos ^ — p sin 9) {X — q cos 9 — p sin 6*) = 0. 



The relation (32) gives the correspondence rule 

q W{q,p) i — )■ ^m9 
Analogously we get the correspondence rule 



d d 



dX 89 



+ X cos 9 w(x, 



pW{q,p)< — -^+Xsm9]wix, 



(33) 



(34) 



The correspondence rules (29), (30), (33), (34) give the possibility to transform the Moyal equation (26) 
into the evolution equation for the optical tomogram w{X, 9). We get the new result in explicit form 

d 



dt 



w{X,9,t) 



+ 2 



Im U < sin 9 



d_ 
09 



dX 



-1 



+ X cos 6* + z 



w{X,9,t) 
sin^ d 



2 dX 



w{X,9,t). 



(35) 



This quantum evolution equation for the probability distribution w{X, 9,t) is compatible with the known 
evolution equation for the symplectic tomogram M(X, /i, z/, t) which can be obtained from the Moyal 
equation in view of the correspondence rules 



q W{q,p) i — > 



pW{q,p) i — )■ 



d \ 


-'A 


dX ) 




d ^ 




dX I 





M{X,fi,u), 
M{X,fi,u), 



d_ 
dq 

d_ 
dp 



d 



W{q,p) M i2—M{X,iJ,iy), 

Wiq,p) ^ zyAM(X,/x,z/). 
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(36) 



These rules are easily obtained by considering Eq.(|9]) where delta function is presented by the Fourier 
integral, i.e. 



M(X,/i,i.) 



47r2 



W{q,p) e^^'^^-^^'-'P^dk dq dp. 



For example the application to both sides of Eq.(37) the operator ^ (^) ^ "we have 



d_f d 
djjt, 



(37) 



(38) 



and equality gives first correspondence rule in (36). Analogously we get other correspondence rules in 



(36). Applying the correspondence rules to Moyal evolution equation we get the evolution equation for 
the symplectic tomogram M{X,fi, u), i.e. 

-M{X,fi,u,t) = fi—M{X,fi,u,t) 



+ 2 





' d ' 


-1 


Im t/j- 






dX 





d iu d 

+ 



M(X,/x,z/,t). 



(39) 



9/i 2 dX 

Using connection of symplectic and optical tomograms ^ one can derive (35) from (39). The 



equation (35) from (39) can be transformed into equations for characteristic functions (16) and (17). 



This can be done by substitutions 



dX 



d_ 
dX 



in (39) in the case of the function xm(-2,/^, J^) and 



_d_ 
dX 



_d_ 
dX 



-1 



X i 



d_ 



Thus we get the evolution equation 



dxM{z,n,i^,t) d 



dt 



U 



\ d vz 
iz djji 2 



Ul- 



Id vz 
iz 9/i 2 



(40) 



for the function Xm{z, /i, z^, t). For the optical tomographic characteristic function we obtain the evolution 
equation 



dt 



d 



sin 29 [1 + 



1] d 
2d^ 



U 



-u 



09 

sin 6* 9 . „ (9 sin 6* 

IT] o9 o9 2 

sin 9 d d sin 9 

IT] o9 o9 2 



(41) 
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4 Multidimensional case 



The results of previous sections can be generalized to multidimensional case. Thus, for multimode 
Hamiltonians 



n „2 

(7=1 



(42) 



and for optical and symplectic multidimensional tomograms 



(X, e,t) = j W{q, p, t) JJ exp - Qa cos 9^ - P^^-^^ | 

/n 
W{q,p, t) Y[ exp{ik„{X^ - fi„q„ - u^p^)} 
a=l 



sin^<,^^ ] (Tk d"g d>_ 
(T k d"g d> 



(43) 



(44) 



(27r;i)2" ' 

where u^j and are frequency and mass dimensional constants respectively for a— th degree of freedom 
{ucr are choosing from convenience considerations for particular Hamiltonian) , multidimensional analogs 



of the equations (|35|) and (|39|) can be written in the forms: 
d 



dt 



w{x,e,t) 



d 1 

cos^ Oa— sin 26 a + X„ sin 6'o- 



de, 2 



d 



+ 



h 



Im U < sin 9^ 



_d_ 



d 



1 -1 



X„ cos 0„ + i 



dX, 

h sin 6 a d 
2maUa- dXa- 



w{x,e,t) 



w{X,9,t)- 



(45) 



d_ 
dt 



M{X,fI,iJ,t) 



CT=1 



M{X,fl,iJ,t) 



+ 



h 



Im U 



d 



dX^ 



d iurrh d 
+ 



dfXa- 2 dX^ 



M{X,fI,i7,t). 



(46) 



We reconstructed in these formulas dimensional constants including the Planck constant. Note, that the 
parameters fi^ are dimensionless, while parameters have dimension of [m]~^[t]. 



For completeness we present also energy level equations for multimode Hamiltonians (42), when U{q) 
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is independent of time. The equations have the form: 



cos^ 9„ 



(7=1 



d 



-2 



2 



' d ' 




_dX,_ 





cos^ + sin 26*0- 



del 

d_ 



+ 1 



^sin^^, 
2 



1 (9^ 

- COS t^, 



+ 



Re U I sin ^, 



d 



d 



E Me{X,II,iJ) 



d9, 

±1- 



dX, 



+ x^cos e„ + I- -— - 

zm^LOa oXa 



2m„ 



d 



dX, 



d 



d 



+ 



Re U 



d 



dX„ 



d ^ iu^h d 



(9/io- 2 dX„ 



dX, 
Me{X,i1,i7). 



"dxi 

WE{X,e); 
ME{X,pi,U) 



we{X, 



(47) 



(48) 



The stationarity conditions of the distributions we{X, 9) and Me{X, fl, z7) are obtained from the equa- 



tions (45) and (46), under conditions dtWE = and dtME = respectively. Thus each tomogram 
of stationary state of quantum system obey two equations simultaneously: energy level equation and 
stationarity condition. 



5 Classical Liouville equation in tomographic form 

Classical Liouville equation in phase space for the potential U{q) has the form 

dfjq, P,t) ^ Pa df{q, p, t) dU {q, t) df{q, p,t) ^ ^ 

dt ^ m„ dq„ ^ dq^ dp^ 

where n - is a number of degrees of freedom. If we introduce optical and symplectic tomograms of the 
distribution function / as follows: 

w,i{X,9,t) = j /(g,p,t)f[exp|zA;. |^X. _ cos^. - p^^^^ | d"g d> d"fc; (50) 



Mc/(X, /I, z7, t) = j^-y^ j f{q, p, t) Y\ exp {iK {X^ - ^^q^ - i^aPa)} d"/c d"g d>, (51) 



the inverse transforms are given by the Fourier integrals: 

+00 



f{q,p,t)= [ / "^ro Zn ^ciiX, 1 1) n ^ exp [iK (x^ - q, cos 9^ - pj-^) ] ; (52) 
J J (27r)2" 1^ V m^^^J } 
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1 /■ - " 

J M^i{X, II, u, t) Yl exp {iK (^a - Aia^a - i^aPa)} d"X d"/i d"i/. 



(53) 



As the relations (50) and (51) for f{q,p) are similar to the relations (43) and (44) for the Wigner 



function W{q,p) we can use multimode analogs of the correspondence rules (29), (30), (33), (34) and 



(36) to obtain evolution equations for the tomograms Wd and Wd from the Liouville equation (42). After 



calculations we get: 
d 



dt 



WdiX,9,t) 



d 



1 



cos^ 9^— - sin 29^ {1 + X, 



+ 



d 



.(7 = 1 



09 „ 



U <qa sin 9o 



d 



dX, 



d 



d 



dX, 



+ X(j cos 9a 



Wd{X,9,t) 

sin 9a d 



rUaUa dXa 



Wd{X,9,t); 



(54) 



-Md{X,fl,u,t) = fl—Md{X,fl,u,t) 



+ 



d 



(7 = 1 



d 



dXa 



d 



d 



d^a f " dXa 



Md{X,fI,u,t). 



(55) 



Note, that these equations are the limit cases of the corresponding equations (45) and (46) when Planck 
constant is taken to be zero. Other classical kinetic equations can be presented in the analogous tomo- 
graphic form. 



6 Tomograms of time dependent photon-added coherent states 
of parametric oscillator 



As an example let us find tomograms of studied in ^2^ time dependent photon-added coherent states 
\a,m,t) of one mode parametric oscillator with the Hamiltonian 

The state \a,m,t) is defined as follows: 

|a,m,t)= U{t)\a,m,0)= U{t)\a,m) = {m\Lm{-\a\'^)y^^^ U{t)a+"'\a), (56) 

where Lm{z) = Lm{z) is the Laguerre polynomial |2T1[22], \a) is the initial coherent state, U{t) is the 
unitary evolution operator 

U{t) U+{t) = i, f/(0) = i. (57) 
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As shown in [20] the expression for the state |a,m,t) in the coordinate representation has the form 
(g|a,m,t) = (m!L„(-|an)"^/^ - ) {q\a,t), 



where Hm{z) is the Hermite polynomial 



exp 



{q\a,t) is the time-dependent coherent state 



2e e 2e 2 
which was found for the first time in [23], c-number function e{t) satisfies the equation 

e{t) + n\t) e{t) = 0, 

with the initial conditions e{0) = 1, £(0) = i, which means that the Wronskian is 

ee* — e*e = —2i. 

The state \a,m,t) is pure, thus we can find tomogram of it from the formular 
M(X,/i,z/,t) = -^ 

After some calculations we get 

(m!L^(-|a|2))"^ 



exp ( - ) {q\a,m,t) dq 



(58) 



(59) 



(60) 



(61) 



Mam.{X,IJ., Z/, t) 



\a\ 



H„ 



'Xe + iV2: 



av 



|e|(/i£: + ve) 



\e^{liE + ve) 
2e^ j \e'^{^e* + ui*) 



1/2- 



exp 



+ 



V2aX a^e* 



iva 



2|/ie + z/e:|2 fie + ue 2e e{jj,e + ue) 



(62) 



The substitutions /i = cos6' and u = sinO to (62) gives us the optical tomogram Wam{X,6,t). 



7 Conclusion 



To summarize, we point out the main results of this work. 

We obtained the evolution equations for optical tomogram of state of quantum and classical systems 
both one dimensional and multidimensional cases. We shown the correspondence between the evolution 
equations for optical and symplectic tomograms. The importance of the quantum tomographic equations 
written namely for optical tomogram is connected with the fact that namely these tomograms are 
measured in experiments to get the photon states and to study their characteristics. The evolution 
equation found in the work provide the possibility of monitoring the system quantum states in the 
process of the time evolution. The generalization of the results to relativistic like in [18], [19] kinetic 
equations will be given in further publications. 
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